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FOREWORD 

This  document  describes  analytical  and  computational  studies  of  three- 
dimensional  incompressible  laminar  and  turbulent  wall  jets  in  small  cross 
flows.  This  effort  was  performed  during  the  period  June  30,  1976,  to 
June  30,  1977,  and  was  sponsored  by  the  Naval  Air  Development  Center  under 
Contract  No.  N62269-76-C-0382 . 

The  technical  monitor  for  this  study  was  Dr.  K.  A.  Green. 
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Modern  naval  aircraft  can  reduce  strike  force  vulnerability  by  distributing 
these  vehicles  over  a larger  number  of  ships  within  the  fleet.  One  way  of 
achieving  this  allocation  is  through  the  attainment  of  vertical  lift-off  capa- 
bility. A technique  used  to  provide  vertical  lift  without  oversizing  the 
engine  in  the  cruise  mode  is  the  use  of  thrust  augmenting  ejectors.  With  these 
devices,  engine  thrust  can  be  enhanced  during  vertical  takeoffs  and  landings. 
Obviously,  it  is  desired  to  achieve  the  highest  thrust  augmentation  ratio  (<p' 
as  possible.  Various  design  concepts  have  been  advanced  toward  obtaining  this 
goal.  In  the  Navy/Rockwell  International  XFV-12A,  for  example,  an  ejector 
system  composed  of  a centerbody  and  two  Coanda  wall  jets  is  currently  under 
development.  A central  feature  of  the  flow  fields  produced  by  this  system  is 
three-dimensionality.  This  has  been  particularly  evident  in  subscale  flow 
visualization  on  the  Coanda  surfaces.  It  is  believed  that  these  flow  processes 
may  be  important  toward  (|)  maximization.  One  way  of  understanding  this  rela- 
tionship is  through  theoretical  modeling  which  can  provide  a means  of  reducing 
the  high  cost  of  powered  lift  testing.  Unfortunately,  existing  methodology 
has  been  limited  in  the  past  to  two-dimensional  flows  for  the  analysis  of  wall 
jets  and  complete  ejector  systems. 

In  Ref.  1,  a semi-analytical  solution  for  a wall  jet  over  a flat  plate  is 
considered.  Both  the  cases  of  laminar  and  turbulent  flow  are  treated. 

Similarity  solutions  are  studied  for  the  laminar  case  in  which  the  flux  of 
exterior  momentum  flux  is  an  invariant.  For  the  flat  plate  case,  the  existence 
of  this  constant  does  not  depend  on  similarity.  With  regard  to  two-dimensional 


1 


j 


J 


Rockwell  International 

Science  Center 

SC5079.5FR 


laminar  jets  over  a curved  wall  treated  in  Ref.  2,  similarity  is  necessary 
to  obtain  the  corresponding  invariant.  Two-dimensional  turbulent  wall  jets 
were  also  considered  by  Giles,  et  al.,^  who  studied  self-preserving  behavior 
for  logarithmic  spiral  profiles.  Various  workers  have  studied  turbulent 
processes  experimentally  in  two-  and  three-dimensional  wall  jet  flows.  This 
effort  is  exemplified  by  Refs.  4-7.  Coflowing  jets  which  in  contrast  to  the 
submerged  case  have  the  jet  embedded  in  an  external  invlscid  field  are  of 
great  practical  interest.  Kruka  and  Esklnazi  have  investigated  deviations 
from  similitude  in  such  flows  as  well  as  merging  of  the  mixing  and  wall  layers. 

Three-dimensional  turbulent  processes  have  been  studied  in  connection 
with  downstream  behavior  of  non-circular  jets  over  flat  plates  and  are 
exemplified  by  Refs.  9 and  10.  These  investigations  have  relevance  to  the 
prediction  of  ejector  three-dimensional  mixing  described  in  Ref.  11. 

The  mathematical  prediction  of  these  flows  presents  formidable  problems. 
Only  the  simplest  geometries,  e.g.,  flat  plate  or  special  wall  shapes  such 
as  the  logarithmic  spiral,  lead  to  an  ordinary  differential  equation  for  a 
similarity  solution.  For  turbulent  flows,  with  realistic  eddy  viscosity 
models,  partial  differential  equations  govern  the  flow  field.  Modern  finite 
difference  methods  offer  promise  of  handling  these  cases.  In  particular, 
Dvorak  in  Ref.  12  treats  two-dimensional  wall  jets  over  boundaries  of  large 
curvature.  Computational  modeling  of  three-dimensional  generalizations  of 
these  flows  has  up  till  now  been  unexplored  to  the  best  of  our  knowledge. 

This  class  of  flows  occurs  in  connection  with  taper  and  sweep  effects  on  lift 
augmenters  and  upper-surface-blown  wings. 
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To  shed  light  on  the  associated  flow  patterns,  a study,  "Three-Dimensional 
Flow  of  a Wall  Jet,"  was  initiated  by  the  Naval  Air  Development  Center  to 
investigate  wall  jet  flows  which  exemplify  typical  features  of  complex  pro- 
pulsive lift  applications.  The  purpose  of  this  study  has  been  to  apply  modem 
computational  methods  to  the  treatment  of  three-dimensional  wall  jets.  The 
following  three  basic  tasks  were  performed: 

Task  1:  Formulate  a model  to  describe  a 3-D  wall  jet  in  the  small 
cross  flow  approximation. 

Task  2:  Develop  a numerical  method  and  computer  code  to  treat  a 3-D 
wall  jet. 

Task  3:  Parametric  studies  using  computer  code. 

In  Task  3,  the  streamwise  developments  of  shear  stresses,  sideslip  angles, 
streamwise,  and  spanwise  velocity  profiles  have  been  studied. 

This  report  will  summarize  the  basic  results  for  all  three  tasks. 
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2.0  FORMULATION  OF  THE  PROBLEM 

2.1  Description  of  Physical  System  and  Assumptions 

The  configuration  shown  in  Fig.  la  has  formed  the  basis  of  this 
investigation.  Depicted  is  a section  of  a three-dimensional  wing  OPCEFO 
which  has  a wall  jet  over  its  surface  ADEF  generated  by  the  efflux  from 
the  slot  ABCD.  An  intrinsic  coordinate  system  (x,y,z)  is  arranged  so  that 
the  slot  ABCD  is  embedded  in  the  surface  x = 0,  and  the  wing  is  the  surface 
y * 0.  Surfaces  x = constant  are  normal  to  the  wing  and  orthogonal  to 
y * constant  as  shown  in  Fig.  lb.  For  simplicity,  a cylindrical  arrangement 
is  shown  with  the  z direction  parallel  to  generators  of  the  cylinder.  How- 
ever, the  formulation  to  be  discussed  can  be  applied  to  more  complicated 
three-dimensional  shapes. 

2.2  Incompressible  Navier-Stokes  Equations 

To  serve  as  a framework  for  subsequent  developments,  the  incompressible 
Navier-Stokes  equations  are  considered  in  this  section. 

Denoting  an  arc  element  ds,  and  the  orthogonal  curvilinear  coordinate 
system  given  in  Fig.  la,  and  the  metric  coefficients  h^ , i “ 1,2,3,  ds  is 
given  by 

2 2 2 2 2 2 2 
ds  - h^dx  + h^dy  + h^dz 

If  u,  V,  and  w are,  respectively,  the  velocity  coordinates  in  the  x,  y and 
z directions,  then  if  q - (u,v,w),  p ■ pressure,  p ■ density. 
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Fig.  la  Geometry  of  wall  jet  configuration 
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Fig.  lb  Intrinsic  coordinate  system 
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U)  ■ vorticlty  ■ curl  q,  then  the  equations  of  motion  for  a laminar  flow* 
with  constant  kinematic  viscosity  v are: 

Continuity 


div  q ■ 0 


Momentum 


qxiu  ■ grad|^  + - V div  grad  q 


On  taking  components,  these  equations  become 


(2.1) 


(2.2) 


Continuity 


(h^h^u)^  + (^3hj^v)y  + (h^h^w)^  - 0 


X Momentum 


2v  2, 

V h-  /■u  \ w ii_ 

2 uu  v(h-u)  3 

X • X t 1 y I w ^ • V 

hihj  hj  hjhj  hjhj  - 


(2.3a) 


(2.3b) 


= 7^.  . 1 

*2^3  u\ 

+ 3 1 

3a\  3 

3x  \ 

hj^  3x  / 

9y 

1 ^2  3y/  3z 

.Turbulent  flows  will  be  considered  in  Section  2.6 

"coordinate  variable  subscripts  indicate  partial  differentiation  with 
respect  to  these  variables. 
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V Momentum 


2 w „ u u 

- 5^  S/  ^ ^ ■ "A  % 


vAv  - 


12  “2  3 


hzP 


(2.3c) 


z Momentum 


(h.w)  + ~ ti.P 

h h 1 h-  h.h,  3 x i y n3"2  3 

3 1z  3 31 


(2.3d) 


2.3  Small  Cross  Flow  Approximation 

Assming  that  «.  3/3g  « 1.  Eq=-  «-3) 

Continuity 

(h,h.u)  + (hjl-^v)  . 0 


X Momentum 


v(h,u) 

4.  ^ 2. 

■ W ^ h,  ^^2 


7 


Rockwell  International 

Scianc*  Center 

SC5079.5FR 


0^ 


Rockwell  International 

Science  Center 

SC5079.5FR 


We  consider  a wall  jet  limit  involving  the  jet  exit  height  H becoming  small 
in  comparison  to  the  wall  radius  of  curvature  at  a fixed  downstream  station 
X.  The  jet  height  5 is  0(H)  as  H 0.*  In  addition,  y,  the  normal  coordinate 
to  the  surface  is  0(6)  in  the  limit.  Furthermore,  we  assume  that  metric 
coefficients  h^  and  (i  ■ 1,2,3)  defined  below  are  0(1).  in  this  limit, 
the  approximate  orders  of  magnitude  of  the  various  terms  are  shown  above  the 
equations  tabulated  below: 


Continuit5 


(h3U)^  + 0.3^)^  . 0 


(2.4a) 


(For  both  terms  to  balance,  v therefore  ~ 6u) 


X Momentum 


uu 

^ (h  u) 

1 y 


^ + vAu 


(2.4b) 


(where  conclusions  from  the  continuity  equation  have  been  used  in  the  ordering) 


Considering  two  arbitrary  functions  f(x)  and  g(x),  f « 0(g)  as  x “*■  x implies 
that  |f/g|  < k as  X x where  k is  independent  of  x.  The  statement  f * 0(g) 
implies  that  f/g  "*■  0 as°x  *♦  0. 
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y Momentian 


2;?  ^ 2 2 
U 0 ou  u 


UV  ^ i. 

X . Ku 

— — + w — r — 

1>1  1 hj 


v6u6 


- -2:  + vA^ 
P 


-2 


(2.4c) 


z Momentum 


u 


uw  uw 


5iiw6 


-1 


^2  2 
0 U 


KjU  + + -j 


UW  v(h,w) 
+ 


X . iZi  4.  K 


vw6 


^ + vAw 
^3 


-2 


(2.4d) 


where 


K,  • h,  /h,h,  , -K_  H h,  /h,h 


‘3  13 

X 


1 ' 13 
z 


= h,  /h. 


’ 1 + Ky  " “1  ' “1 

7 


-K3  S /h^h^ 
z 


2.5  Finite  Momentum  Limit  for  Finite  Curvature  Walls,  (<  » 0(1))  — 
Submerged  Wall  Jets 

To  further  simplify  the  foregoing  equations,  we  consider  the  limit  in 

y 

which  pu  £ is  fixed  as  6 0.  Here,  u ■ 0(n) , where  U is  defined  as  a mean 

—1/2 

jet  exit  velocity.  Accordingly,  u ■ 0(5  ).  If  < ■ 0(1)  as  5 “*■  0,  h^^  = 1, 

allowing  various  temns  to  be  eliminated  from  the  foregoing  equations.  The 
ground  rules  for  this  process  are  that  at  least  the  frictional  term  in  the 


More  complex  forms  of  the  equations  arise  for  <6  ■ 0(1)  but  will  not  be 
considered  in  this  report. 
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streamwlse  momentxim  equation  is  retained,  and  a nontrivial  y momentum  is 


tions  become,  no 


z Momentum' 


hi IL 


. 
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Here,  the  most  interesting  case  has  been  retained  so  that  the  K^u  term 


balances  other  terms  in  the  0(uw)  z momentum  equation  since  it  is  of  order 
u^j^  which  itself  is  assumed  to  be  0(uw) . For  submerged  wall  jets  with 


0,  as  y in  contrast  to  curved  wall  boundary  layers,  the  pressure 


gradient  term  p /ph.  is  negligible  in  (2.4b),  since  from  (2.4c), 


p - 0(<u  5)  * 0(1).  A similar  result  is  obtained  in  the  finite  mass  limit 


Pu6  “ fixed  as  5 “*■  0.  Only  for  boundary  layers  or  subregions  of  coflowing 


wall  jet  flows  with  11m  u ^0  jets  Implying  p = 0(u  ) in  (2.4b)  can  the 

y~*co  ^ ^ 


streamwise  pressure  gradient  become  important.  For  the  finite  momentum 


2/3 

limit,  inclusion  of  the  friction  term  in  (2.4b)  implies  6 ~ v as  v “*■  0. 


This  order  of  magnitude  has  been  tacitly  assumed  in  the  omission  of  the  higher 


order  term  vv  in  (2.5c). 

yy 


The  rationale  for  the  6 scaling  with  v and  the  disappearance  of  the  p 


term  from  the  x momentum  equation  for  submerged  jets  can  be  more  fully  under- 
stood from  three-dimensional  generalizations  of  asymptotic  developments  to 
be  discussed  shortly  in  connection  with  two-dimensional  flows.  Prior  to 
this,  we  note  that  for  finite  mass  with  pu6  fixed,  5 ~ v as  V “*■  0.  As  will  be 
indicated,  other  "distinguished  limits"  are  possible  in  which  internal  struc- 
tures such  as  the  wall  layer,  potential  core,  and  mixing  layers  can  be 
abstracted. 

We  conclude  this  section  by  noting  that  the  foregoing  approximate  forms 
of  the  equations  of  motion  could  be  obtained  from  a formal  asymptotic  expansion 
procedure  which  will  be  illustrated  for  two-dimensional  curved  wall-jets.  It 
is  well  known  that  these  flows  can  be  divided  into  a transitional  region  near 
the  jet  exit  consisting  of  a mixing  layer,  inviscid  constant  velocity  potential 
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I 

f : 

I 


i 


I 


li 


I 

i 


core,  and  a boundary  layer  in  the  vicinity  of  the  wall.  The  potential  core 

is  eaten  up  by  the  boundary  and  mixing  layers.  Turbulent  diffusion  results 

in  a merger  of  these  layers  at  a downstream  location.  In  what  follows,  we 

consider  the  flow  in  the  fully  merged  zone.  As  in  boundary  layers,  two 

different  representations  can  be  used  to  describe  the  flow  structure.  An 

"inner"  representation  is  appropriate  to  the  viscous  jet  layer  near  the  wall, 

and  an  "outer"  expansion  describes  the  external  inviscld  flow.  Another  option 

is  to  develop  a uniformly  valid  asymptotic  representation  using  an  optimal 

13  14 

set  of  coordinates  developed  by  Kaplun.  ’ 

Denoting  the  mean  velocity  at  the  exit  by  U = Q/pH,  where  H is  the  exit 

height,  and  Q is  the  exit  mass  flow,  the  exit  momentum  is  QU.  Accordingly, 

the  nondimensional  form  of  Eqs.  (2.4)  in  two  dimensions  can  be  obtained  by 

normalizing  all  velocities  with  respect  to  U,  the  pressure  difference  from 

2 

ambient  with  respect  to  pU  , and  all  lengths  with  respect  to  L a mean  radius 

fc 

of  curvature.  The  resulting  dimensionless  equations  of  motion  are  similar 
in  form  to  (2.4)  except  with  suitable  dimensionless  redefinitions  of  the 
p and  the  v coefficients  replaced  by  R where  R * Reynolds  number  based  on 
L - UL/v. 

We  now  consider  appropriate  asymptotic  representations  for  the  inner 
viscous  layer.  Introducing  a small  parameter  e which  is  the  reciprocal  of 
the  Reynolds  number  R,  we  envision  a sequence  of  flows  observed  at  a fixed 
X station  in  which  the  normalized  wall  height,  H,  is  allowed  to  become 
vanishingly  small  as  e “*■  0.  If  5(x;e)  is  the  characteristic  jet  height,  6 

* I 

Note  that  other  normalizing  lengths  are  possible  such  as  the  viscous  length 
v/U  or  the  jet  height  H.  The  velocities  can  also  be  referred  to  a free 
stream  velocity  provided  the  latter  is  not  zero.  The  selection  mode 
here  is  advantageous  for  the  arguments  that  follovr. 
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will  scale  like  H as  H 0.  To  keep  the  fine  structure  of  the  jet  layer  in 
view  as  H ^ 0,  we  "blow  up"  the  y scale  by  a factor  o(e),  where  the  functional 
foirm  o(t)  is  to  be  determined.  Here,  a ~ 5.  To  formalize  this,  we  assert 
that  the  y dependence  is  really  a dependence  on  the  strained  variable  y = y/o. 
The  most  general  form  of  the  inner  expansion  leading  to  the  non-dimensional, 
laminar  two-dimensional  specialization  of  Eqs.  (2.5)  is 


u(x,y;e) 

-2  ~ -1 
= ea  u^(x,y)  + to  + . . . 

(2.6a) 

v(x,y;e) 

= £0  ^v  (x,y)  + ev,  + ... 

0 1 

(2.6b) 

p(x,y;e) 

2-4  2-3 

= £ o p (x)  + e a p (x,y)  + . . . 
o 1 

(2.6c) 

for  an  "inner  limit,"  x,y  fixed  as  e i 0.  The  "gauge  function,"  a is 
determined  by  matching  this  solution  with  the  outer  inviscid  flow. 

It  should  be  recognized  that  for  the  flat  plate  boundary  layer,  since 
there  is  no  characteristic  length  in  the  strearawise  direction,  the  appropriate 
representations  are  coordinate  expansions  for  large  x rather  than  for  small 
values  of  the  parameter  £ of  (2.6).  Another  viewpoint,  see,  for  example. 

Van  Dyke,^^  is  to  introduce  a fictitious  normalizing  length  in  the  streamwise 
direction  which  cancels  out  in  the  analysis. 

For  wall  jets,  several  "distinguished  limits"  are  relevant  for  the  co- 
flow ratio  Ci)  = U/U^,  where  is  the  freestream  velocity  in  the  outer  flow. 
These  cases  are  as  follows : 

(i)  0)  “*•  0 

(ii)  0)  fixed 
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(lii)  (J  ^ 

(iv)  (0  =•  “* 

as  e "*•  0.  Case  (i)  is  not  of  interest  for  propulsive  lift  applications. 

Note  further  that  Case  (ii)  subsumes  Case  (iv)  which  corresponds  to  a sub- 
merged jet.  If  Case  (ii)  is  assumed,  then  the  assertion  that  u = ^1),  uni- 
formly in  0 < X ^ is  plausible  based  on  normalization  of  this  streamwise 

velocity  component  with  respect  to  U and  matching  considerations.  Accordingly, 
*"2 

ea  = 1 in  (2.6a)  implying  that  a = v^.  This  scaling  is  also  appropriate  to 
conventional  boundary  layer  flows.  Substitution  of  (2.6)  into  the  exact 
equations  and  retaining  terms  of  dominant  order  will  give  the  non-dimensional 
analog  of  (2 .5a)-(2.5c) , for  the  approximate  quantities  in  (2.6),  with  an 
additional  pressure  gradient  term  in  the  axial  momentum  equation  due  to  the 
coflow  effect.  These  equations  are: 

Continuity 


u + V =0 
o o~ 

X y 


X Momentum 


u u + V u 
o o o o~ 
X y 


-p'  (x)  + u 
o o-~ 

yy 


Momentum 


■="0  ’ n- 
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The  longitudinal  gradient  is  determined  as  in  conventional  boundary 

layers  by  matching  with  the  outer  flow  streamwise  pressure  gradient  which  is 
determined  from  Bernoulli’s  equation.  Note  that  the  y pressure  gradient 
balancing  centrifugal  force  across  the  streamlines  arises  from  the  second 
order  term  p^^  in  the  pressure  expansion  (2.6c). 

The  representation  of  the  outer  flow  field  is  obtained  from  other 
asymptotic  expansions  of  the  flow  variables.  The  appropriate  outer  variable 


i: 

i: 


normal  to  the  body  surface  is  y and  the  expansions  are: 


u(x,y;e)  = U^(x,y)  + t/e  U^(x,y)  + ... 

v(x,y;e)  = V (x,y)  + V (x,y)  + ... 
o 1 

p(x,y;e)  = P^(x,y)  + /i  Pj^(x,y)  + ... 


for  x,y  fixed  as  e 0 ("outer  limit"). 

On  substitution  of  these  expansions  into  the  exact  equations  and 
retaining  the  dominant  terms , the  following  equations  are  obtained  for  the 
first  order  quantities: 


Continuity 


U + (h.V  ) = 0 

o 1 o y 

X 


X Momentum 


U U + h,V  U + KU  V - -P 
oo  loo  oo  o 

X y y 
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-3 


i 


U V + h.V  V 
o o loo 
X y 


To  detetmlne  the  longitudinal  pressure  gradient  In  the  Inner  equation  and 
a(e),  Bernoulli's  equation 


-2 

0) 


and  a matching  procedure  Is  used  in  which  the  Inner  and  outer  solutions  are 

written  in  a representation  appropriate  to  an  intermediate  "overlap"  domain 

between  inner  and  outer  regions  in  which  the  solutions  have  common  validity. 

For  this  purpose,  the  intermediate  limit,  y^  fixed  as  e 0,  is  used  in  which 

y^  “ y/Tl(G)  and  the  order  of  ri(e)  is  between  /e  and  unity.  The  inner  and 

outer  expansions  are  written  in  terms  of  y^  and  are  equated  to  various  orders, 

* 

yielding  conditions  on  the  unknown  quantities . 

From  the  Bernoulli  equation  and  this  procedure,  the  following  boundary 
conditions  are  obtained 


f ii'  V (x,0)  “0  on 

U' 

\ ■ P^(x,0)  - Pq(x)  - 0)  ^ - U^(x,0)/2 


Van  Dyke  in  Ref.  14  uses  Lagerstrom's  restricted  matching  principle  to  obtain 
I similar  results  for  boundary  layers  without  the  intermediate  variable 

f formalism  applied  in  this  section.  Cole  in  Ref.  17  has  applied  the  inter- 

I mediate  variable  matching  method  for  a wide  class  of  singular  perturbation 

'i  problems  and  has  derived  formulations  similar  to  those  described  here  for 

/ boundary  layers. 
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U (x.O) 
o 


u (x,®°)  = u (x) 
o e 


Pj^Cx.y) 


Vj^Cx.O) 


(x)  - 2u  (x)U,  (x,0)  as  y -►  <» 
e el 


5*'(x)  , 5* 


= f [u  - u 

Jo  ° « 


The  solution  procedure  is  to  solve  the  outer  equations  with  the  first 

of  the  above  boundary  conditions.  The  quantity  u is  subsequently  used  with 

e 

to  solve  the  inner  problem  with  an  initial  condition  of  the  form 


Uo(0,y)  = g(y) 


where  g is  a prescribed  function.  In  this  respect,  and  the  turbulence  models 
employed,  the  wall  jet  problem  differs  from  the  boundary  layer  formulation. 
The  latter  derives  its  initial  conditions  from  matching  with  the  outer  flow, 
whereas  for  wall  jets,  these  are  specified  independently. 

The  remaining  boundary  conditions  comprise  the  no-slip  conditions, 
u^(x,0)  = v^(x,0)  = 0 for  all  i,  and  the  outer  boundary  conditions  for  p^ 
appearing  in  the  inner  y momentum  equation. 

Note  that  the  quantity  Uj^(x,0)  must  be  obtained  from  the  solution  of 
the  second  order  outer  problem  with  V^(x,0)  expressed  in  terms  of  the  slope 
of  the  displacement  thickness  6*'(x) . Higher  approximations  are  obtained 
using  a similar  iteration  procedure  relevant  to  this  weak  viscous  Interaction 
problem. 


18 


Rockwell  International 

Sci«nce  Center 

SC5079.5FR 


Aside  from  the  differences  noted,  the  foregoing  problem  strongly 
resembles  a boundary  layer  formulation,  for  Cases  (ii)  and  (iv) . The  latter 
case  is  obtained  from  the  former  by  letting  u^  = = 0 for  all  i. 

The  longitudinal  pressure  gradient  term  in  the  inner  x momentum  equation  is 
thereby  eliminated. 

If  the  velocities  are  normalized  with  respect  to  for  Case  (iii) , 
scalings  for  the  inner  variables  are  obtained  which  correspond  to  those 
derived  in  the  dimensional  formulation  given  in  Section  2.5.  Formally,  the 
inner  equations  become  in  this  case 


-1/3  / ~x  ^ 1/3  ^ 

u ■ e u (x,y)  + e u,  + . . . 
o 1 

V * e^^^v  + ev,  + . . . 

o 1 

-2/3  ~ 2/3 

p « e ' ^ ^ + ... 

^ 2/3 

for  x,y  = y/e  fixed  as  e "*■  0.  This  will  yield  identical  inner  equations 
to  dominant  order  as  for  Cases  (ii)  and  (iv).  However,  it  is  anticipated 
that  details  of  the  matching  will  be  different.  As  a check,  Rickley's 
similarity  solution  for  a free  jet  with  transverse  momentum  flux  invariant 
along  the  jet  exhibits  the  same  e scaling  shown  in  the  dominant  terms  of 
the  foregoing  expansions. 

It  is  noteworthy  that  the  submerged  jet  of  Case  (iv)  is  degenerate  with 
respect  to  (ill) . This  is  plausible  since  normalizations  of  the  latter  are 
non-existent  for  U„  • 0. 
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2.6  Turbulence  Assumptions — Eddy  Viscosity  Models 

In  previous  studies  conducted  at  the  Science  Center,  a number  of  turbu- 
lence models  were  investigated.  Because  of  the  orientation  of  this  investiga- 
tion to  algorithm  development,  a detailed  study  of  the  adequacy  of  these  models 
was  not  attempted.  However,  it  should  be  noted  that  the  numerical  algorithm 
to  be  described  in  subsequent  sections  is  general  enough  to  assimilate  the 
various  turbulent  models  which  for  the  purpose  of  the  present  investigation 
have  been  restricted  to  eddy  viscosity  simulations.  For  purposes  of 
discussion  of  the  niimerical  algorithm  and  the  results,  the  turbulent  frame- 
work corresponding  to  Eqs.  (2.5)  differs  inthe  respect  that  the  terms 
and  vw^^  in  the  laminar  formulation  are  replaced,  respectively,  by  their 
eddy  viscosity  counterparts  ((v+Cj^)u^)y  and 

A prototyplc  model  selected  to  illustrate  the  application  of  a typical 
eddy  viscosity  simulation  is: 


(0.435y)‘ 


3y 


(0.125yj^)^ 


3u 


3y 


. y > y 


* 


(2.7a) 


(2.7b) 


where  y^  is  determined  by 


u(yj^)  - 0.01 

Uy(yi)  < 0 


It  should  be  noted  that  this  model  provides  coupling  between  the  spanwise 
flow  w and  the  streamwise  field  u not  occurring  in  the  weak  cross  flow 
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laminar  formulation.  From  a computational  viewpoint,  the  coupling  was 
suppressed  to  achieve  an  efficient  algorithm.  In  our  procedure,  the  finite 
difference  approximation  used  for  (2.7)  is  such  that  the  discretized 
momentum  equations  are  effectively  decoupled.  This  was  achieved  by  evaluating 
3v/3y  at  the  previous  streamwlse  station  Instead  of  evaluating  the  average 
between  the  present  and  last  computed  streamwlse  station. 

Other  turbulence  models  have  been  proposed  for  two-dimensional  wall  jets 
in  which  the  wall  curvature  affects  the  entrainment  and  eddy  viscosity  simu- 


lation. These  can  be  accommodated  by  our  computational  procedure. 
2.7  Boundary  and  Initial  Conditions 


The  boundary  conditions  to  be  employed  are  the  no-sllp  conditions  at 
the  wall  and  as3nnpcotlc  conditions  relevant  to  an  "outer"  flow  field  external 
to  the  jet.  Thus,  on  the  wall  y ■ 0, 


u(x,0,z)  * v(x,0,z)  ■ w(x,0,2)  « 0 


(2.8a) 


At  the  edge  of  the  jet,  y ■ «» 


u(x,'*»,z)  - u^(x,z) 


(2.8b) 


v(x,“,z)  - V (x,z) 
e 


(2.8c) 


w(x,«',z)  • w^(x,z) 


(2.8d) 


p(x,“»,z)  - P.(x,z) 


(2.8e) 
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The  initial  profile  must  satisfy  the  compatibility  conditions,  i.e.,  a solu- 
tion of  (2.5)  or  its  turbulent  counterpart  evaluated  at  x * 0,  subject  to 
the  appropriate  specialization  of  (2.8).  It  should  be  noted  that  in  an 
incompressible  context,  the  quantity  p^  can  be  determined  from  Bernoulli's 
theorem  providing  the  outer  flow  is  invlscld. 

2.8  Formulation  in  Glauert  Variables 


To  minimize  sharp  gradients  and  smooth  the  computational  problem,  the 
governing  equations  of  motion  are  rewritten  in  a new  set  of  independent  and 
dependent  variables.  The  Glauert  wall-jet  transformations  given  in  Ref.  1 
are  used  to  change  the  Independent  variables  (x,y)  to  (s,n)*: 


ds  ■ hj^dx 


(2.9a) 


1 -3/4 

n - 4 3 y 


(2.9b) 


A new  dependent  variable  f(s,ri)  is  introduced  such  that 


1/4 

^ * s ' h2f(s,n)  , 


(2.9c) 


where  ’■p  is  the  stream  function  satisfying  the  continuity  equation  with 


The  quantities  (u,v,x,y ,1^;)  are  dimensionless,  being  obtained  from  the 
corresponding  dimensional  variables  (u,V,R,V,-J;)  by  writing  3 - Uu, 

V ■ Uv,  X • vx/U,  y - vy/U,  ij)  ■ v^'p/U. 
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Furthermore,  let  p be  the  reduced  pressure  given  by 


p * As^'^S/P 


(2.9d) 


Using  these  transformations,  the  equations  of  motion  (2.5)  simplify  to: 
Streamwise  Momentum  Equation 


ff  + 2r  - 4s(f  f - f^f  ) 
nn  n u us  s nu 


(2.10a) 


Vertical  Momentum  Eqtiation 


p - <f 


(2.10b) 


Spanwise  Momentum  Equation 


* ^2>“n)n  * {^*^-3r)%*  ‘'4  V * 


K-f  • 4s (fw  - f w ) 
2 n ns  s n 


(2.10c) 


In  addition,  the  boundary  conditions  are  given  by 


Boundary  Conditions  at  the  Wall,  n ~ 0 


f(s,0)  - f (s,0)  • wCs.O)  - 0 


s > 0 


(2.11a) 
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Boiindary  Condition  at  Jet  Edge 


R(s) 


w(s,oe)  ■ W(s) 


(2.11b) 


p(s,oo)  - P(s) 


where  R and  W are  arbitrary  fxinctlons  of  s obtained  from  the  external  flow, 
and  P can  be  obtained  from  Bernoulli's  theorem.  Consistent  with  the  small 
cross  flow  approximation,  the  dependence  on  z is  absent.  It  is  tacitly 
assumed  in  the  streamwlse  momentum  equation  that  P^(s)  « 1,  otherwise,  the 


equivalent  of  the  “P^/p  term  should  be  added  to  the  right-hand  side  of  the 


streamwlse  momentimi  equation  in  accord  with  a three-dimensional  qualitative 
extension  of  the  matching  procedures  elucidated  in  Section  2.5. 


I 
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3.0  NUMERICAL  METHODS 


3.1  The  Box  Scheme 

To  solve  the  wall-jet  equations  in  Section  2.7,  an  implicit  finite 

18 

difference  method  (the  Box  Scheme)  developed  by  H.  B.  Keller  is  used. 
The  differential  equations  are  written  as  a first  order  system  in  terms 
of  relabeled  dependent  variables  u(s,n),  v(s,il).  t(s,ri): 


(3.1a) 


U ■ V 

n 


(3.1b) 


Wfi  - t 


(3.1c) 


(d  + 


-f  ^ it)  * *"<“"3  ■ 's''> 


(3. Id) 


4K, uw  - K„u  + 4s (uw 
12  s 


(3.1e) 


p ■ <u 


(3. If) 


(In  passing,  we  note  that  the  right-hand  side  of  the  above  system  (3.1)  does 


not  Involve  terms  that  are  derivatives  of  u.) 

I In  I I d 

Now  consider  any  family  of  meshes  2 they 

satisfy  the  following 
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I 
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s *0 
o 


la  -s  ,+k  . n«l,2,...,N 

^ n n-1  n ’ 


, edge  of  jet 


The  following  notations  are  used: 


Vl/2  ■ 2 Vl> 


1j-l/2  ■ 2 '"j  "j-l' 


n-1/2  1 


■j  2 ■ “n-" 

‘j-1/2  - 


Tlj)  + Z(S 


Tlj)  + Z(S 


Vr"  ■ • 


i-i')  ’ 


s_,n  - n 


C.1  - 1 + Cl  , 


0-2  - 1 + Cj 


1 Wz  3»/j. 


n-1/2 
1/2 


n-1/2 


k 
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We  now  derive  the  difference  equations  approximating  system  (3.1).  From 
Fig.  2,  consider  box  Eqiiations  (3.1a-c)  are  approximated  by 

centering  about  segment  1*2^3  streamwise  station 

at  i^lch  the  solution  vector  (f ,u,v,w,t,p)  is  to  be  computed): 


j 1-1  n 

hj  “ '^j-1/2 


n n 

_J izl  = n 

h^  ''j-1/2 


n n 

w.  - w.  , 

-1 Izi  » t°  , 

hj  ''j-1/2 


(3.3a) 


(3.3b) 


(3.3c) 


Next,  Eqs.  (3.1d-f)  are  approximated  by  centering  about  * 


the  middle  of  the  box  P.P.P.P, : 

12  3^ 


(a,v)“-(a,v)“ 


V']  '''^l''M-l  „ „ 2.n  / c"-!  n ^n  n-1  \ 


-n  n \ n-1 

j-1/2  ■ ^j-l/2'^j-l/2/  ^1 


(3.3d) 


*^6(''j-l/2“j-l/2  'j-l/2^J-l/2^  * '2 


L /r,  / n n ,n  -n  \ 

■*■  ^^6'“j-l/2''j-l/2  " ^j-l/2'^j-l/2' 


(3.3e) 
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~n  -n 


/ .11-1/2/  2.n  , n-1 

^'^^j-l/2^“  ^j-1/2  *3 
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(3.3f) 


where  g”,^,  and  g^"^  (the  dependent  variables  in  g^,  are 

evaluated  only  on  the  previous  streanwise  station  ate  given  by 

.n-l  ^n-l 


n-1 


j 

. / / 2. n-1  , ,n-l  n-1  \ 

+ 4Pg(-(u  ^j-l/2^j-l/2' 


(3.3d^ 


g 


n-1 


.n-1 


n-1/2/  .n-1 


(1  + (ft)j_j^/2  ■ ^^^^j-l/2^'^''^j-l/2 


,„,.n-l/2/  2, n-1  , / n-1  , n .n-1  ..^n-1  /^n-l  ,n  .\ 

Vl/2^^2^J-l/2^'^  ^j-1/2  ^^6("''j-l/2^'^j-l/2'^j-l/2^  ^j-l/2^^j-l/2"^j-l/2^) 


(3.30*) 


n-1 


Si 


/~.n-l  ~n-l 
(P)^  P 


‘‘j 


+ ^'^''j-1/2^'"  ^j-1/2 


Equations  (3.3a-3.3f)  together  with  (3.3d-f^)  are  to  be  applied  to  all 
n-points,  j » 1,2,...,J. 

The  boundary  conditions  to  be  applied  at  s » s are: 

n 


wall 


n 


w„ 


(3.4a) 

(3.4b) 
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(3.4c) 


(3.4d) 


3. 2 Solution  of  the  Difference  Equations 

Assuming  solution  Is  known  at  s = s ,,  l.e.,  u?~^. 

n-1’  1 J j ’ 

~n“l\  £ r.  j X 

Wj  » » Pj  ) for  0 ^ j ^ J,  we  now  want  to  evaluate  the  solution  at 

s “ s^.  We  apply  Eqs.  (3.3)  for  j = 1,...,J.  Together  with  the  boundary 

conditions  (3.4),  this  yields  6*(J+1)  equations  for  the  6*(J+1)  unknowns 

-n  n n n n n.  . „ , 

^j’''j’''j’"j’’'j’^j^’  ^ ” 0,1,2,...,J. 

For  turbulent  wall  jets,  the  streamwlse  and  spanwise  momentum  equations 
are  coupled  through  the  eddy  viscosity  of  Eqs.  (2.7).  To  handle  this  compu- 
tationally, the  streamwlse  momentum  equation  is  solved  using  the  9w/3y  asso- 
ciated with  the  previous  s step  as  indicated  in  Section  2.6.  This  effectively 
decouples  the  streamwlse  momentum  equation  from  the  spanwise  momentum  equation 

at  any  station  s • s , reducing  the  computational  time  and  storage  requirement. 

n 

With  the  above  approximation,  the  solution  algorithm  is  then  given  by: 


(1)  Solve  for  (fj,u°,Vj)  , j - 0,1,2,...,J 
(11)  Use  (1)  to  solve  for  (Wj.t”)  , j - 0,1,2,...,J  (3.5) 

(ill)  Use  (1)  to  compute  p”  , j « 0 ,1 ,2 , . . . , J 


, .L  IjHlWUlPPiW Jll  JJW.I) 


<!» 
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To  treat  (i),  we  must  solve  a system  of  3*(J+1)  nonlinear  equations.  (The 
equations  are  (3. 3a,b,d)j_j^  and  (3.4a,c).)  Then,  assuming  (f,u,v)  is 
successfully  computed  at  s ■ s^,  (il)  and  (lii)  involves  only  systems  of 
linear  equations.  Thus,  the  major  bulk  of  computational  time  is  in  (i) . 

We  note  that  the  difference  equations  (3.3)  for  the  variables 
(f°,Uj ,v^,Wj ,t°,p^)  for  j ” 0,1,2,...,J  and  s ■ s^  can  be  viewed  as  the 
solution  to  two-point  boundary  value  problems  of  systems  of  linear  or  non- 
linear ordinary  differential  equations,  with  the  independent  variable  being 
r|.  Thus  (i)  now  can  be  viewed  as  solution  to: 


df 

— • u 

dn 

du 

dn  " "" 


d(Ct.v) 

dn 


-d  + P^)fv  - 2u^  + 4Pg(u^  - fv)  + g^(n) 


(3.6a) 

(3.6b) 

(3.6c) 


with  boundary  conditions 


f(0)  - u(0)  - 0 
u(oo)  ■ constant 


(3.6d) 

(3.6e) 


we  have  deliberately  suppressed  the  dependence  of  s in  P^,  P^,  g^,  and  the 
constant  in  (3.6e).  However,  they  do  change  as  we  march  downstream. 


31 


^<r4ggg  *«tgyr j j 


vlv  Rockwell  International 

^ Sci«nc«  Canter 

SC5079.5FR 

The  theory  of  numerical  solution  to  two-point  boundary  value  problems 
for  ordinary  differential  equations  can  be  found  in  Refs.  19  and  20.  We 
shall  only  outline  the  procedure  here. 

_ ..n  n n.  J 

The  nonlinear  system  of  eqxiations  for  the  unknown  U = ’^j' j«o 

to  be  solved  by  Newton's  method.  Specifically,  we  define  ^(U),  (suppressing 
the  s-dependence  in  (f,u,v)  again). 


"o-^i 


V"o 

• 2 

V +v 

0 1 

‘^i 

2 

a-v,-a 

V 

1 1 ( 

2-^+  a 

Uj  - constant 


(3.7) 
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Let  the  initial  iterate  be  the  solution  at  the  previous  streamwise 

I * 

station  s , . Then,  Newton  s method  gives 
n~i 


|i  ^ , V > 1 


3U 


(3.8a) 


u(v)  . 


, V > 1 


(3.8b) 


Method  is  said  to  have  converged  at  the  Kth  iteration  when 


^prescribed 


error  tolerance 


(3.9) 


The  Jacobian  matrix  3$/3D  in  (3.8a)  has  a very  nice  structure,  a consequence 
of  the  centered-Euler  method  in  approximating  (3.6) 


0 (010)-' 


*When  depends  on  , 2.  * K-1,K,  as  most  all  eddy  viscosity  models  do,  then 
we  do  riot  have  Newton's  method  strictly  speaking,  because  we  avoid  terms 
[ (3aj^/3u)iSu]v  and  [ (3ot^/3v)6v]v. 


i 


I. 
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where  are  (3x3)  matrices,  K » are  given  by: 


with 


^2  - i W-1/2  ^"6'^-l/2} 

®3  " ■ 2 {^^■‘■^l^^K-1/2  ^^6^K-l/2} 


(3.11a) 


(3.11b) 


(3.12) 


The  first  two  rows  In  (3.10)  are  contributions  from  boundary  conditions  at 
the  wall,  with  the  last  row  from  the  jet  edge.  3$/9U  can  be  further 
partitioned  into  a block  tridiagonal  matrix  [B.A. C. ],  where 
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i > 1 


(3.13a) 


i > 1 


(3.13b) 


i > 1 


(3.13c) 


Before  we  go  Co  the  next  section  to  describe  an  algorithm  for  solving  such 

a matrix  system,  we  want  to  comment  on  the  solution  procedure  for  (11)  and 

(111)  of  (3.5).  As  remarked  earlier,  since  (f,u,v)  are  now  known  at  s ■ s , 

n 

the  equations  for  ^ “ 0,1,2,...,J  can  be  viewed  as  the  difference 

approximation  to  the  linear  two-point  boundary  value  problem  In  one  Independent 
variable  n of  the  form 


^ - A(n)z  + 3(Ti) 

(3.14a) 

B^(0)  - 0 

(3.14b) 

* constant 

(3.14c) 
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T 

with  Z = (w,t)  , A has  coefficients  (f,u,v)  and  g (O.g^),  the  system  of 
linear  equations  can  again  be  partitioned  into  the  form 


where  A.  is  a block  tridiagonal  matrix,  Z^  and  b^^  are  given  by 

0 
0 

gzChi/z) 

0 

*2^^j-l/2 


®2^'^J-l/2 
0 

The  computation  of  (ill)  is  simply  the  centered-Euler  integration  of 

-OO 

p(n^)  - p(®®)  - / <U  (T)dT 
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Ajj,  Cjj  are  matrices  of  order  n,  K ■ 1 J. 

A can  be  decomposed  Into  the  form 


A - JCu  5 [3^  I 0]  * [0  c^] 


(3.18) 


where 


“l"*! 


(3.19a) 


«iVi  • “l 


, 1 * 


(3.19b) 


ct^  • “ ^1^1— 1 * i“2,3,...»J 


Here  matrices  in  turn  are  decomposed  into  the  form: 


cti  - 


(3.20) 


where  p^,  are  permutation  matrices  for  row-and-column  pivoting.  and 

are  lower  and  upper  triangular  matrices.  It  is  important  to  have  an 

accurate  LU-factorization  of  because  they  are  used  In  solving  both 

and  x^.  Here  we  use  a mixed  pivoting  strategy.  During  the  Kth  stage  of 

(K) 

Gaussian  elimination,  the  pivot  aj^j^  is  chosen  to  satisfy 


4 > k 


(3.21) 


38 


it'  ' ' 


•^^BBgWigWBBPIIliIjillipiil  II.  JAU.-.  j- 
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This  mixed  pivoting  Is  much  better  than  the  partial  column  pivoting 


or  partial  row  pivoting.  If  e.  Is  the  round-off  error  using  the  mixed 

m 

pivoting  strategy,  and  If  Is  the  round-off  error  using  either  partial 
column  or  partial  row  pivoting  strategy,  then  It  can  be  easily  shown  that 


le,  II  >lle,  II 

Ip  Im 


(3.22) 


The  solution  of  from  Eq.  (3.19b)  can  be  easily  carried  out.  Assuming 
there  are  p rows  of  with  at  least  one  non-zero  element  on  that  row,  then 
the  solution  of  8^^  corresponds  to  Inverting  the  following 


~T  _T 

’"i-iZk  - ® 


^1-1% 


K 


1,2,. ..,p 


(3.23) 


where  lower  and  upper  triangular  matrices  of  size  n.  We 

further  note  the  zero  structure  (as  in  (3.13a))  is  preserved  under  such  a 
decomposition  scheme.  This  avoids  unnecessary  storage  space  requirements. 

Now  assume  (3.19)  has  been  performed,  then  to  solve  x,  we  merely  have 
to  solve 


(3.24) 
C3.25) 

T T 

where  z = . . . ,Zj)  , Zj^  ■ {z^^z^,, ..  ^z^  can  be  obtained: 


JCz  ” b 
Ux  - z 
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h " 


~K  * ~K  ■ ®K~K-1 


K - 2,3,. .. ,J 


(3.26a) 


(3.26b) 


and  (3.25)  will  give  the  solution  vector  x 


(3.27a) 


'^l-l^l-l  “ ^£-1  "*  ^ “ J,J-1,...,2 


(3.27b) 


3.4  Starting  Procedure 

The  starting  procedure  Is  to  employ  a suitable  discretization  of  the 
Initial  conditions  obtained  by  the  method  described  In  Section  2.7.  Glauert’s 
similarity  solution  derived  In  Ref.  1 was  Implemented  using  the  s >■  0 
soe'.iallzatlon  of  Eqs.  (3.3)  and  (3.4)  for  the  results  given  in  this  report. 
However,  with  the  compatibility  restrictions  given  in  Section  2.7,  more 
general  initial  conditions  could  be  accommodated  including  those  derived  from 
experimental  data. 


1 
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In  this  section,  results  will  be  indicated  typifying  calculations  using 
the  computational  model.  Because  of  its  interest  on  the  XFV-12A  augmenter, 
a logarithmic  wing  spiral  contour,  shown  schematically  in  Fig.  3,  will  be 
discussed.  In  contrast  to  previously  published  solutions  exemplified  by 
Ref.  3,  this  discussion  will  deal  with  non-similar  flows  due  to  the  nature 
of  the  assumed  turbulence  model.  In  Ref.  3,  the  logarithmic  spiral  shape 
with  certain  assumptions  on  the  scaling  of  .jet  thickness  with  downstream 
distance  gave  rise  to  similitude  and  an  analytic  solution  for  the  flow.  The 
non-similar  framework  considered  here  makes  such  a solution  unlikely  and 
numerical  methods  must  be  used.  In  the  notation  of  the  figure,  the  equation 
of  the  spiral  contour  is 


9/K 


(4.1) 


where  s is  the  running  arc  length,  9 is  the  local  inclination  of  the  surface, 
where  K and  Sq  are  constants.  For  K > 0,  a convex  contour  is  obtained,  and 
with  K 0,  concavity  is  implied.  Equation  (4.1)  can  be  represented 
parametrically  as 


(4.1a) 


(4.2b) 
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In  Figs.  7 and  8,  the  streamwise  development  of  the  u and  w profiles  is 
shown.  Although  both  profiles  resemble  each  other,  the  momentum  in  the  cross 
flow  increases,  in  contrast  to  the  decay  exhibited  by  u.  This  trend  is  also 
indicated  in  Fig.  9 for  w„,_  and  is  due  to  the  source-like  manner  in  which 
the  sidewash  is  produced. 

To  further  assess  the  influence  of  turbulent  coupling  of  the  sidewash 


field  on  the  mainstream  flow,  the  effect  of  w on  typical  velocity  profile  is 


LOG  SPIRALS 


spirals  for  K = .05  and  K 


development  of  various  wall  jet  flow  quantities 


= 0,-5,  f(0,«>)  » 1 
Submerged  Case,  » 0 

p = Physical  Pressure 


\ 


/-Similarity  Solution  Slope 


4Vs” 


Fig.  7 Streanwlse  development  of  u profiles  for  log 
spiral-submerged  wall  let  with  scan  flow 


Rockwell  International 

Science  Center 

SC5079.5FR 


[}J^r. 


shown  in  Fig.  10,  where  the  peak  region  is  magnified  to  show  the  very  small 
effect  of  the  crossflow. 

To  illustrate  the  resemblence  of  velocity  profile  of  coflowing  wall 
jets  and  conventional  boundary  layers,  a point  made  in  Section  2.5,  calcula- 
tions were  performed  using  the  typical  model  of  Eqs.  (2.5)  with 
p ■ w • Kj  ” f(0,0)  ■ f^(0,0)  ■ 0,  and  uCx,®*)  » 1.  Results  for  the  stream- 
wise  development  of  the  reduced  velocity  profile  and  shear  stress  on  a flat 
plate  with  fCO,**)  * 4 are  shown  in  Figs.  11  and  12.  To  indicate  the 
potentialities  of  the  existing  code,  the  streamwlse  velocity  profile  develop- 
ment with  downstream  distance  is  shown  in  Fig.  13  for  a logarithmic  spiral 
contour  with  cross  flow.  Here,  K2  ■ -10,  uCx,®®)  - f(0,o“)  - 1,  and  the  external 


pressure  gradient  was  neglected  in  the  calculations.  With  this  assumption 
the  qualitative  downstream  behavior  resembles  that  of  a flat  plate. 


I 

f 


Fig.  10  Effect  of  spanwise  flow  on  normalized  u profile 


La  51 


I 3.0  CONCLUSIONS  AND  RECOMMENDATIONS 

i A computational  model  using  Keller's  box  scheme  has  been  developed 

to  treat  incompressible  turbulent  wall  jets  in  a small  cross  flow  approxima- 
tion. The  coiiq)uter  code  can  handle  sidewash  w injected  as  a source  term 
in  the  spanwise  momentum  equation.  The  effect  of  the  span  flow  on  the 
streamwise  flow  is  due  to  the  eddy  viscosity  coupling  between  the  u and  w 
fields.  For  this  type  of  spanwise  flow  generation,  the  coupling  appears 
extremely  weak,  reducing  the  peak  streamwise  component  and  causing  growth 
; of  w momentum  in  the  downstream  direction. 

In  subsequent  effort,  different  modes  of  sidewash  addition  will  be 
Investigated,  i.e.,  through  the  boundary  and  initial  conditions.  The  results 
of  this  are  indicative  of  flow  conditions  for  wall  jets  on  configurations 
with  slight  taper  or  sweep.  For  handling  more  realistic  situations,  the 
foregoing  analysis  will  be  extended  to  finite  cross  flow. 


Rockwell  International 

Science  Center 

SC5079.5FR 


6 . 0 REFERENCES 

Hi.  B.  Glauert,  "The  Wall  Jet,"  J.  FI.  Mech.  1,  625-643  (1956). 

2 

I.  J.  Wygnowski  and  F.  H.  Champagne,  "The  Laminar  Wall  Jet  Over  a Curved 
Surface,"  J.  FI.  Mech.  Ft.  3,  459-465  (1968). 

J.  A.  Giles,  A.  P.  Noyes,  and  R.  A.  Sawyer,  "Turbulent  Wall  Jets  on  Logarithmic 
Spiral  Surfaces,"  Aero  Quarterly  XVII,  202-215  (1965). 

4 

R.  J.  Goldstein  and  D.  J.  Wilson,  "Turbulent  Wall  Jets  with  Cylindrical 
Streamwise  Surface  Curvature,"  ASME  Trans.,  J.  Fluids  Eng.,  550-557  (1976). 

^H.  P.  Irwin  and  P.  A.  Smith,  "Prediction  of  the  Effect  of  Streamline 
Curvature  on  Turbulence,"  Phys.  of  Fluids  (6),  624-630  (1975). 

^R.  M.  C.  So,  "A  Turbulence  Velocity  Scale  for  Curved  Shear  Flows," 

J.  FI.  Mech.  70,  Pt.  1,  37-57  (1975). 

^I.  P.  Castro  and  P.  Bradshaw,  "The  Turbulence  Structure  of  a Highly  Curved 
Mixing  Layer,"  J.  FI.  Mech.  7^,  Pt.  2,  265-304  (1976). 

Q 

V.  Kruka  and  S.  Eskinazi,  "The  Wall  Jet  in  a Moving  Stream,"  J.  FI.  Mech. 

20,  Pt.  4,  555-579  (1964). 

Q 

N.  V.  Chandrasekhara  Swamy  and  P.  Bandyospadhyay , "Mean  and  Turbulence 
Characteristics  of  Three-Dimensional  Wall  Jets,"  J.  FI.  Mech.  7^,  Pt.  3, 

541-562  (1975). 

^^P.  M.  Sforza  and  G.  Herbst,  "A  Study  of  Three-Dimensional  Incompressible 
Turbulent  Wall  Jets,"  AIAA  J.  8,  276-283  (1970). 

^^A.  D.  De  Joode  and  S.  V.  Pantankar,  "Prediction  of  Three-Dimensional 
Turbulent  Mixing  in  an  Ejector,"  AIAA  Paper  77-706  (1977). 

12 

F.  A.  Dvorak,  "Calculation  of  Turbulent  Boundary  Layers  and  Wall  Jets  Over 
Curved  Surfaces,"  AIAA  J.  3^,  517-524  (1973). 


56 


Rockwell  International 

Science  Center 

SC5079.5FR 


Kaplun,  Z.  Angew  Math  u Mech. , 111-135  (1954). 

A.  Lagerstrom,  "Laminar  Flow,"  B,  from  Laminar  Flows  and  Transition  to 
Turbulence,  Princeton  Series  on  High  Speed  Aerodynamics  and  Jet  Propulsion. 
Vol.  IV  (Princeton  University  Press,  1962). 

Van  Dyke,  "Higher  Approximations  in  Boundary  Layer  Theory,"  Part  1, 
General  Analysis,  J.  FI.  Mech.  161-177  (1962). 

^^M.  Van  Dyke,  Perturbation  Methods  in  Fluid  Mechanics  (Academic  Press, 

New  York,  1964),  p.  124. 

D.  Cole,  Perturbation  Methods  in  Applied  Mathematics  (Blaisdell, 

Waltham,  Mass,  1968). 

B.  Keller,  "A  New  Difference  Scheme  for  Parabolic  Problems,"  Numerical 
Solutions  of  Partial  Differential  Equations,  Vol.  II  (1971). 

B.  Keller,  "Accurate  Difference  Methods  for  Nonlinear  Two-Point  Boundary 
Value  Problems,"  SIAM  J.  Num.  Anal.  11,  305-320  (1974). 

B.  Keller,  "Accurate  Difference  Methods  for  Linear  Ordinary  Differential 
Systems  Subject  to  Linear  Constraints,"  SIAM  J.  Num.  Anal.  6_,  8-30  (1969). 


57 


Rockwell  International 

Science  Center 


SC5079.5FR 


APPENDIX  A;  CODE  LISTING  AND  SAMPLE  INPUT  AND  OUTPUT 
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Sample  Computer  Output 
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APPENDIX  B:  3-D  WALL- JET  SMALL  CROSS  FLOW  PROGRAM-RUNNING  INSTRUCTIONS 

The  following  deck  setup  will  indicate  control  cards  and  data 
(a  name  list  card  INPUTS)  needed  to  nm  the  following  example: 

(A)  Small  cross  flow  CASEW^.TRUE. 

(B)  u ■ 0 at  outer  edge  of  jet,  C4  ■ 0. 

(C)  Parameter  in  logarithmic  spiral  K » 1/3,  Cl  ■ 1./3. 

(D)  The  Induced  magnitude  of  induced  cross-flow  * 1,  C7  » 

(E)  Meshes,  initial  profile  + streamwlse  output  stations  are 
be  provided  by  program. 

(F)  Streamwlse  station  to  be  solved  to  XB  « 1. 

(G)  Initial  mass  flux  > 1. 

1.  Deck  Setup 
Job  Card 
FTN4. 

LGO. 

7/8/9 
Source  program 
7/8/9 

-$  INPUTS  CASEW-.TRUE.,  C1-. 33333333,  C4-0.,  C7— l.$ 

t 

Second  column 
6 / 7 / 8 / 9 

For  ocher  desired  cases,  see  definitions  of  the  various  variables  and 
their  options  In  the  listing.  The  solutions  are  to  be  printed  on  paper 


cards 
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In  sequence  of  streamvlse  stations.  Thus,  for  any  s location,  there 
are  eight  columns  of  outputs: 

Column  1:  Index  of  vertical  net  points 

Column  2:  y — vertical  net  point  values,  from  y ■ 0 (1  In  column  1, 
at  the  wall)  to  y « YB  (the  last  value  In  column  1,  outer 
edge  of  jet) 

The  next  six  columns  have  the  same  convention  as  column  2 : 

Column  3:  F — Glauert  similarity  variable 

Column  4:  DF — the  partial  derivative  3/3 (eta)  F 

Column  5:  DDF — the  partial  derivative  3/3 (eta)  DF 

Column  6:  P — the  reduced  pressure  P 

Column  7:  W — cross-flow  velocity 

Column  8:  DW — the  partial  derivative  3/3 (eta)  W 

To  find  out  the  computed  values  of  (F,  DF,  DDF,  P,  W,  DW)  at  y ••  0.5, 
say,  we  need  to  look  at  the  horizontal  line  with  the  y-value  on  column 
2 to  match  with  y * 0.5  (provided  y ■ 0.5  is  a net  point),  then  the 
third  to  eighth  columns  will  give  the  function  values  of  F,  DF,  DDF, 

P,  W,  DW  at  y ■ 0.5. 
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